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Abstract 

We study the [i —>- ey decay in the ^-invariant next-to-minimal supersymmetric (SUSY) Standard 
Model (NMSSM) with superheavy right-handed neutrinos. We assume that the soft SUSY breaking 
parameters are generated at the GUT scale, not universally as in the minimal supergravity scenario 
but in such a way that those soft parameters which are specific to the NMSSM can differ from 
the soft parameters which involve only the MSSM fields while keeping the universality at the GUT 
scale within the soft parameters for the MSSM and right-handed neutrino fields. We call this type 
of boundary conditions “semi-constrained”. In this model, the lepton-flavor-violating off-diagonal 
elements of the slepton mass matrix are induced by radiative corrections from the neutrino Yukawa 
couplings, just like as in the MSSM extended with the right-handed neutrinos, and these off-diagonal 
elements induce sizable rates of fi —> ey depending on the parameter space. Since this model has 
more free parameters than the MSSM, the parameter region favored from the Higgs boson mass can 
slightly differ from that in the MSSM. We show that there is a parameter region in which the /i —> ey 
decay can be observable in the near future even if the SUSY mass scale is about 4 TeV. 


1 Introduction 


It is now clear that the lepton flavor number is not a conserved quantity because of experimental ob¬ 
servations of neutrino oscillations [T]. In the minimal extensions of the Standard Model (SM) with the 
Majorana neutrino mass terms, the branching ratios for charged lepton-flavor violating (LFV) processes 
are extremely small since they are suppressed by at least a factor of m 2 /ni \ v , which makes it very diffi¬ 
cult for near-future experiments to detect LFV signals. On the other hand, in more general extensions 
of the SM, which are motivated by several reasons, it is known that sizable LFV rates are predicted 
depending on parameter region. If LFV processes are discovered, it directly means an indirect signa¬ 
ture of physics beyond the SM (BSM). Recently, the MEG experiment reported a new upper limit of 
Br(/i — > ey) < 5.7 x 1CU 13 [2]. This already gives a strong constraint on models beyond the SM, and 
hence it is very important to keep updating these upper bounds on the LFV processes. 

Supersymmetry (SUSY) is still a promising candidate for physics beyond the SM [3]. Lots of effort has 
been devoted to the discovery of SUSY at the LHC, but only in vain so far. The most studied model of 
SUSY is the minimal SUSY SM (MSSM). Even in the framework of the MSSM, there are some unsolved 
problems such as the /r problem. Next-to-the MSSM (NMSSM) is an extension of the MSSM with a 
SM-singlet Higgs chiral superfield S. The NMSSM could give a hint to solve the fi problem since in this 
model the /x term is induced by the vacuum expectation value (VEV) of the scalar component S of S. In 
this sense the NMSSM is a natural extension of the MSSM. 

One of the difficulties in the MSSM is the Higgs boson mass. In the MSSM, the tree-level lightest 
Higgs boson mass is bounded from above as, 


m h.MSSM 


< M§ cos 2 2/3 


( 1 . 1 ) 
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and has to rely on large radiative corrections to reproduce the observed Higgs boson mass of 126 GeV pQ. 
The main contribution to the radiative corrections comes from the top Yukawa coupling 0HS], and to 
maximize this effect one needs a top-squark mass much larger than the top-quark mass. In the NMSSM, 
the lightest Higgs boson mass reads [7]: 

m h, nmssm ~ M z c °s 2 2/3 + sin 2 2/3 + Am 2 h 1Loop , (1.2) 

where v ~ 174 GeV. As is seen from this equation, the contribution from the new parameter A, which is 
the coupling among the new singlet S and the MSSM Higgs doublets H u and Hd , makes the tree-level 
Higgs boson mass larger, in particular for small tan/3. We have to note that the mixings between S and 
the MSSM Higgs doublets can make a negative contribution to the lightest Higgs boson mass, and the 
NMSSM does not always predict a larger Higgs boson mass. We will discuss this issue in details later in 
this paper. 

There are more than one-hundred free parameters in the MSSM. Usually, we assume an underlying 
scenario for SUSY breaking, and it allows us to reduce the number of free parameters. In this paper we 
assume the minimal supergravity (mSUGRA)-like boundary conditions that the SUSY breaking param¬ 
eters mo, Mi/ 2 , Aq are universal at the GUT scale. The parameters at the SUSY scale are obtained by 
evolving these parameters according to the renormalization group equations (RGE). These mSUGRA-like 
boundary conditions are very effective for avoiding constraints from the SUSY-induced flavor changing 
neutral current (FCNC) processes. This is also true for the charged LFV processes, and in the mSUGRA, 
also known as the constrained MSSM (cMSSM), there are essentially no charged LFV. This is similar in 
the case of the constrained NMSSM. 

The neutrino masses are exactly zero in the framework of the SM, which clearly needs modifications in 
view of the observation of neutrino oscillations. One of the most natural mechanisms to explain the tiny 
neutrino masses is the (type-I) seesaw mechanism |51U0j. which we consider in this paper. The extension 
of the original seesaw mechanism to SUSY models is straightforward. In the MSSM extended with the 
right-handed neutrinos vr, which we call the MSSM + vr model, even if one assumes the mSUGRA-like 
boundary conditions at the GUT scale, off-diagonal elements in the slepton mass matrices are induced 
via radiative corrections from the neutrino Yukawa couplings, which can predict sizable rates for the 
LFV processes like /i —> eq. This mechanism also works in the NMSSM extended with the right-handed 
neutrinos, which we call the NMSSM + vr model and which we consider in this paper, but since there 
are more free parameters than in the case of the MSSM + vr model, the predicted LFV rates can slightly 
differ from those in the MSSM + vr model in the parameter region favored from the Higgs boson mass. 

The contents of this paper are as follows. In Section 2, we introduce the model we work with, and 
in Section 3 we explain the origin of the LFV (off-diagonal) elements of the slepton mass matrices. In 
Section 4, we discuss constraints on the parameters of the model. We introduce the results of numerical 
calculations in Section 5, and in Section 6 we summarize this paper. 


2 Model 

2.1 Z3-invariant NMSSM 

The NMSSM is an extension of the MSSM, and it has an extra Higgs chiral superfield S which is singlet 
under the SM gauge group. In the Z 3 -invariant NMSSM [7], the /i term /tH u ■ Hd in the superpotential of 
the MSSM is replaced by the term A SH U ■ Hd, and the /r-parameter is determined from the singlet VEV 
s as /x e ff = As. Namely, the superpotential of the Z 3 -invariant NMSSM is given as 

Wnmssm = Hmssm| m=0 + XSH U ■ Hd + -acA 3 , (2-1) 

where the dot in the term A SH U ■ Hd represents the S'C/(2)-invariant product of two SU( 2) doublets, and 
the hats on the fields stand for the superfields corresponding to the fieldtQ. We assume that the A-parity 
is conserved, and assign the even A-parity to S. The soft SUSY breaking terms are 

Isoft = Vmssm| ^—R—o + ^A A\SH U ■ Hd + —kA k S 3 + h. c.^j + TOglS) 2 . (2-2) 

^df we are to consider more general NMSSM, not the Z 3 -invariant NMSSM, we will have more terms in the superpotential, 
such as the terms linear and quadratic in S [7], and also the soft SUSY breaking terms associated with them. To reduce 
the number of free parameters, in this paper we only consider the Z 3 -invariant version of the NMSSM+z/r model. 
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Q 

U c 

D c 

L 

E c 

N c 

H u 

H d 

S 

Z 3 charges 

c 2 

1 

1 

w 2 

1 

1 

UJ 

UJ 

UJ 


Table 1: The Z 3 -charge assignment in the Zg-invariant NMSSM extended with the right-handed neutrinos. 
In the table, u = e 27 ”/ 3 . 


In the case of the constrained NMSSM, the gaugino masses, sfermion soft SUSY breaking masses, and 
the A-parameters take the values which are “universal” at the GUT scale, similarly to the case of the 
cMSSM: 


M a 

— Ml/ 2 


(2.3) 

™ 2 h u 

= m H d = 

ml , 

(2.4) 


= m 0 5ij 

{f = Q,U,D,L,E) , 

(2.5) 

(Af)ij 

o 

II 

(.f = U,D,E) , 

(2.6) 


where a (a = 1,..., 3) labels the gauge groups of the SM, and i and j are the indices for generations, 
i,j = 1,... ,3. As for the parameters A \, A K and m| which are specific to the NMSSM, we assume 
that the values of A\ and A K at the GUT scale are not necessarily equal to Aq 7 and that m| at the 
GUT scale can be different from mg. We call the NMSSM with this class of boundary conditions the 
semi-constrained NMSSM. 

2.2 Z 3 -invariant NMSSM extended with right-handed neutrinos 

In this paper we take the simplest extension of the Zg-invariant NMSSM with the right-handed neutrinos, 
in which the (type-I) seesaw mechanism [SttlOj is at work. The superpotential is given by 

W = Wnmssm + {Y N )jiH u ■ LiNj + ±(M N )ijN?N? , (2.7) 

where the Zg-charges are assigned as in TablefflUIl- This charge assignment excludes the term (A v )i 3 SNf- 
N? from the superpotentiaH. 

The neutrino masses in this model is 

K)ij = u 2 sin 2 /3(Y N )ki((M N )~ l )ki(Y n )ij 

= (t^MNs)*fe m ^,fc(^MNs)fcj j (2-8) 

where Umns is the MNS matrix m and m Ut k (fc = 1,..., 3) are the eigenvalues of the left-handed neutrino 
mass matrix . In the standard representation of the PDG, the matrix reads: 

( C12C1S Si 2 Ci 3 

— Sl 2 C 23 — Ci 2 SigS 2 3 e l<5 Ci 2 C 2 3 — Si 2 Si 3 S 2 3 e i<5 
si 2 s 2 g — c 12 s 13 c 23 e lS ci 2 s 23 — Si 2 si 3 c 23 e lS 

x diag(l,e ia2l/2 ,e M3l/2 ) , 

where c-i 3 = cos Oij , Sij = sin Oij. The mixing angles 0,j ( i,j = 1,..., 3, i < j) describe the mixing between 
the mass eigenstates Vj and v 3 , and the factors S,a 2 \,a 3 \ are complex phases, and represent the Dirac 

2 It is possible to derive the (left-handed) neutrino masses via the type-I seesaw mechanism from the Majorana masses 
which emerge from term (A v)ijSN? • TV? after replacing S with its VEV. In this case, since the singlet VEV ( S) is at most 
0(1 — 100 TeV), the Majorana masses must be about the same order, which forces us to assume a very small neutrino 
Yukawa coupling (Yn) in order to explain the tiny neutrino masses. This makes the LFV rates extremely small and hence 
we do not consider this scenario in this paper. We should also note that in some of extended seesaw schemes such as 
the inverse seesaw mechanism m , it is possible to have the neutrino Yukawa couplings of 0(0.01) and the right-handed 
neutrinos of the mass of the order of the EW scale simultaneously. It is known in literature that in the SUSY inverse seesaw 
models the predicted LFV rates can be sizable m 


—i8 



(2.9) 
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phase and the two Majorana phases, respectively. According to the latest data [I] the values of the angles 
are: 


sin 2 2012 = 0.846 ±0.021 


sin 2 2023 = 


0.9991q'q5s (normal mass hierarchy) 
1 . 000 ^o!oi 7 (inverted mass hierarchy) 


sin 2 2013 = 0.093 ±0.008 . 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


The mass-squared differences, which are also important parameters, are: 

Amli = 7.53 ±0.18 (10“ 5 eV 2 ) 

f 2.52 ± 0.07 (10 -3 eV 2 ) (normal mass hierarchy) 

32 2.44 ±0.06 (10~ 3 eV 2 ) (inverted mass hierarchy) . 

In this paper, we assume the normal hierarchy scenario for the neutrino masses, and take the values 


(2.13) 

(2.14) 


and, for the mixing angles, 


TO„.i 

= 1(T 6 eV 



(2.15) 

m„, 2 

= + 

A m 2 21 

= 0.0087 eV 

(2.16) 


= \/ m l,2 + 

A m§ 2 

= 0.050 eV 

(2.17) 

S12 = 

- 0.55, S23 = 

= 0.66, 

•sis = 0.15 . 

(2.18) 

we take 





S = «21 : 

= OL 3i = 

= 0 , 

(2.19) 


for simplicity. Another free parameters are the 3x3 elements of Mjv- Although it is known that the 
structure of this matrix gives an influence to the predicted LFV rates H5HH, in this paper we assume 


(Mjf)ij = M u x Si, 


( 2 . 20 ) 


where M v is a real number. 

Under the assumption of Eq. (12.201) . we can determine the neutrino Yukawa coupling matrix Y/v from 
Eq. (12.81) by using the input data Eqs. (12.101) (12.191) . How exactly to do this is well known in literature 
(see e.g. [ID])) and below is a brief summary of the procedure. First, in the basis where the charged-lepton 
Yukawa coupling matrix is diagonal, Yy can be expressed by using two unitary matrices U and V and a 
diagonal matrix Y® as 

(Y N ) ii =V ik {Yg) kt U ij . (2.21) 


When the right-handed neutrino Majorana mass matrix is diagonal, by a suitable redefinition of the right- 
handed neutrino superfields, we can take the matrix V to be a unit matrix without loss of generality: 

(Y N )ij = {Y$)iU i0 . (2.22) 

where (Yjy )* (i = 1,..., 3) are the diagonal entries of the matrix namely, ('Yjv)» j = (Xn )Aj- By 

substituting this into Eq. (EH), we can identify U as C/mns and determine (Yy )-; to be 


that is, 


(Y£)i = 


I 47 1 , j 
v 2 sin 2 /3 ’ 


( Y N)ij = 



(UMNs)ij ■ 


(2.23) 


(2.24) 


Later in this paper we use this expression to calculate Br(/i —> eq). 
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3 Lepton Flavor Violation 

In this section we discuss charged lepton flavor violation, taking the NMSSM + vr model as an example 
of new physics beyond the SM. 


3.1 11 —> e 7 in the Standard Model with vr 


Within the SM, the neutrinos are strictly massless and lepton flavor number is exactly conserved. The 
experimental observations of neutrino oscillations [I], however, make it clear that we have to extend 
the SM in such a way that it can accommodate the neutrino masses and mixings. One of the simplest 
extensions is to introduce right-handed neutrinos (hr) which are singlet under the SM gauge group, which 
allows us to introduce Dirac mass terms for the neutrinos in the Lagrangian. 

Once we introduce the right-handed neutrinos in the SM, in general, charged lepton flavor number is 
no longer conserved. This is similar to the case in the quark sector, and the mismatch between the gauge 
eigenstates and the mass eigenstates violates the lepton flavor number conservation. The branching ratio 
of fi —^ ey in this model is given by |2TH23l 


Br(/r —> ey) 


a 

32n 


(t^MNs)*i(^MNs)/xi 

i=2,3 


Mil 


2 

v,l 


2 


(3.1) 


The suppression factor (m^ — i)' 2 /M^ v makes the branching ratio extremely small, and it is very 

difficult for near future experiments to detect p —> ey in this model. On the contrary, in the non-minimal 
extensions of the SM such as the (N)MSSM+iyR, sizable LFV rates can be predicted depending on the 
parameter region, and this makes the LFV searches very important as a probe of new physics beyond the 
SM. 


3.2 fi —> ej in NMSSM + ur Model 

In the NMSSM+zzr model, there are two diagrams which give dominant contributions to the k —> l : f j 
decays (where i and j are the generation indices which run from 1 through 3 with i > j). One is the 
diagram with the neutralino and the charged slepton in the loop, and the other is the diagram which 
involves the chargino and the sneutrino. In general, the amplitude T for the U —> lj 7 decay can be written 
as 

T = emi i e a *Uj(p - q ) [ia a/3 q p (A^Pl + A^Pr)] Ui(p) , (3.2) 

where e is the positron charge, e“ is the polarization vector of the photon, Ui and Uj are the spinors 
for the initial- and final-state leptons, respectively. The momenta p and q are the incoming momentum 
of the initial state lepton li and the outgoing momentum of the final state photon, respectively. The 
operators Pl,r stand for the chiral projection operators. The dependence of the amplitude on the models 
is included in the coefficients A% and A In the case of the MSSM+jzr model, the explicit forms of A% 
and A 2 are given, for example, in Refs. [SSJIMIES] • In the case of the NMSSM + vr model, they are 
essentially the same as the MSSM + vr model, except that there are five neutralinos, instead of four, at 
low energies, and we can use the expressions in Refs. [2312HG03 with small modifications. By using the 
formulas mentioned above, the decay branching ratio Br (li —> lj 7 ) can be calculated from the amplitudes 
to be 

2 ^ 

Br(^ —> lj 7 ) = I 2 + l^?| 2 ) > (3-3) 

where Fq is the total decay width of the lepton li. 

A comment on the singlino (i.e., the fermionic component of the singlet Higgs chiral superfield S ) is 


rrder here. The neutralino mass matrix AI^o 

in the NMSSM is [7]: 



/ Mi 

0 

-grvd/VZ 

giv u /V2 

0 \ 


0 

m 2 

g2Vd/ V% 

-g2v u /Vz 

0 

II 

O 

3 



0 

Meff 

-\v u 





0 

-A v d 


\(symm.) 
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where “symm.” means that this matrix is symmetric, A4T 0 — A4^o, anc ^ Meff = v « — ^sin/3 and 

Vd = ucos/3 where /? = arctan ((H u )/(Hd)). The parameters M\ and M 2 are the U(Y)y and SU(2)l 
gaugino masses, respectively, and g\ and g 2 the U{1)y and SU{2)l gauge couplings, respectively. Since 
for our sample parameters discussed in Section 5 the value of A is taken to be as small as 0.1, and since 
the value of the (5,5) component of -Ad^o for our sample parameters is 0(1 — 10)TeV, which is much 
larger than \v u and A Vd, the mixing between the singlino and the other components of the neutralinos 
are suppressed by the small values of the (3,5), (4,5), (5,3) and (5,4) entries of A4^o. In addition, since 
the singlino does not couple to the MSSM matter fields at tree level, the contribution from the singlino 
component to the LFV rate is negligible for our sample parameters. 

In order to have a non-vanishing LFV rate, we must have off-diagonal elements in the slepton mass 
matrices. The mass matrices are given as, 


Mr = 


M 2 ll 

Mrl 


m 2 lr 

M RR 


= 


= m 


L,ij 


-Mj cos 2/3 <5 


where M| L , M^ R , Mj R . M RL are the 3x3 matrices whose (i,j) elements are given as 

= ml.ij+Vd (^eYe).^ +M|cos2/3(-i +sin 2 6» w )d y - , 
( M RR)ij = m E,ij + v 2 d (yE Y E^.. ~ M% cos 2/3 sin 2 O w Sij , 

{MIr).. = v d ([A* E ) i:j + /rtan/3j ( Ye ^ , 

( m rl) ij = ((MIr) 1 ).. . 


(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 
(3.10) 


In this paper, we assume mSUGRA-like boundary conditions, in which all the SUSY breaking parameters 
that have flavor indices do not have flavor mixings at the GUT scale. This means that there are no off- 
diagonal elements in the matrices M? and M|. However, off-diagonal elements in these mass matrices 
are induced by radiative corrections at the energy scale higher than Mjv, which can be seen in the RGE, 


16lr2 ^(™'i)i, = 


1671-2 


NMSSM 

2vt 


+ ( Y N Y N m L)ij + ( m L Y N YN ) i 3 Y ^( Y N^ m N ) Y n) 

+ 2(1 N Y N)ijVn 2 Ru +2(TjyTjv)jj , 


(3.11) 


where t = InQ with Q being the renormalization scale and (Tjv)y (i,j = 1,...,3) are the trilinear 
coupling among the right-handed sneutrino v Rl the left-handed slepton L , and the Higgs field H u , 

Uneutrino trilinear — ( Y N^)jiH u ■ LiUrj . (3.12) 

The RGE above directly means that both M~ and M? have off-diagonal elements at low energies. The 
size of these off-diagonal elements can be roughly estimated as [5011241125] , 

(A ml)ij = --j^ln M ^ T (6 ml + 2A%)(Yj f Y N ) ij , (3.13) 

where i ^ j- As is clear from Eq. (13.131) . the slepton off-diagonal elements in this model comes from 
the neutrino Yukawa couplings, Yyr. The branching ratio can be estimated in terms of the off-diagonal 
elements to be [20] 


Br(Zj -)■ Ijj) 




((m 2 L h) 2 

Ms USY 


(3.14) 


At this moment, the most stringent experimental constraint on the /i —> ey is given by the MEG exper¬ 
iment and the upper limit is 5.7 x 10' 13 m ■ This bound will be further improved by the upgraded 
MEG experiment to ~ 6 x 10 -14 [50], and this makes the experiment very important as a probe of new 
physics beyond the SM. 
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3.3 Other cLFV processes 

In this paper we focus on fi —»• in the later sections, but there are many other charged LFV pro¬ 

cesses m- Here we mention some of them. 

There are two other U — > Ij'y processes, r —> yi'y and r —> e'y. Their current experimental limits are 
Br(r —> e'y) < 3.3 x 10 -8 and Br(r —> yi'y) < 4.4 x 10 -8 [II. In the near future, these limits are expected 
to be improved to the level Br(T —> I'y) < 1.0 x 1CU 9 at Belle-II [21]. Under the assumptions we set out 
at Section 2, the yi e 7 decay is more sensitive to SUSY particles, and hence we focus on yi —> e'y in 
this paper. 

Other important cLFV processes include if —> if if IJ and yi-e conversion in nuclei. As for the former 
process, when the photon mediation diagram is dominant, the branching ratio can be related to that of 
the k —> Ij'y decay as [ 20 ] 


Br {It ^ If If IJ) 


a /16 mi, 

87T \ 3 2 mi } 



Bi(k -y Ij'y) , 


(3.15) 


and hence Br(Z+ if if l ■ ) can be calculated once Br (k —> Ij'y) is obtained. The current experimental 
limit for — > e + e + e~ is Br(/x + —> e + e + e _ ) < 1.0 x 10 ~ 12 [T|, and this is expected to be improved to 
Br(// + —> e + e + e~) < 1.0 x 10 ~ 16 at the Mu3e experiment at PSI [22]. Concerning the yi-e conversion in 
nuclei, there is a simple relation between the conversion rate B^^N) and Br(/i —> e'y) when the photon 
mediation diagram gives the dominant contribution [30], 


Bfj, e (N) = R(Z)Br(yi -> e 7 ) , (3.16) 

where B fle (N) = T(yi~N — ► e~ N)/T{yi~ N —> capture) is the conversion rate normalized to the muon 
capture rate N capture), and R(Z) is a parameter which depends on the atomic number Z of the 
nucleus which captures the muon. The current limits are B^ e { Ti) < 4.3 x 10 -12 , B^ e ( Au) < 7 x 10 -13 [Tj. 
The near future experiments are the COMET experiment at J-PARC [53]; the Mu2e experiment at 
FNAL [32]; and the PRISM/PRIME experiment at J-PARC [35], which are expected to improve the 
bounds to R /(e (Al) ~ 6 x 10 -17 [54], B^ e ( Al) ~ 6 x 10 -17 [55], B^ e (Al) ~ 10 ~ 18 [55] . respectively. 
Since the R{Z) factors for these experiments are R ~ 0.0025 for Al and R ~ 0.0040 for Ti [30;, these 
experiments are expected to go beyond the corresponding limit of the /i —> e'y decay by 1.5 ~ 3 orders of 
magnitude, and this will be very useful to probe broader parameter region of new physics. 


4 Constraints on the Parameters in the Model 

In this section we discuss constraints on the parameters in the NMSSM + vr model. Some of the issues 
below are already discussed in literature 7\. 

Tadpole conditions 

In the NMSSM, there are three tadpole conditions. At tree-level they read: 

v u + Meff + ^ 2v d + 4 ~ V d) S j ~ VdHefiBe ff = 0 , (4.1) 

v d (ff^Hd + Meff + A-U 2 + 4 ( v d — V u)^j — VuldefiBeS = 0 , (4.2) 

s (m 2 s + kA k s + 2 k 2 s 2 + A 2 (n 2 + vf) - 2A nv u Vd) - A v u VdA\ = 0 , (4.3) 

where yi e s = As and B e s = A\ + ns. We can use these relations to determine three parameters from 

other parameters. For example, we can use these relations to determine /j, e ff,B e s and mf from the other 

parameters. Later we will discuss which parameters we use as input. 

Maximal Tree-level Higgs Mass condition 

One of the advantages of the NMSSM over the MSSM is that there is a parameter region in which the 
lightest Higgs boson mass can be made larger than that of the MSSM. As can be seen from Eq. (11.21) . 
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in order for the Higgs boson mass to be larger, it is favorable to have large A and small tan /3. The 
approximate formula Eq. m is obtained by neglecting the mixings between the MSSM Higgses and the 
singlet Higgs in the CP-even Higgs-boson mass matrix, 

/ Mf cos 2 /3 + HeffBes tan/3 (X 2 v 2 - ^Mf) sin 2/3 - A (2 n eS v d - {B eS + ns) v u ) \ 

■Ms,Tree = M Z shj2 P + Meff^eft COt /3 A (2 ^L eS V u - (B eS + KS ) V d ) , 

\ (symm.) ks (A K + 4ks) + XA\ ) 

(4.4) 


where v u = v sin /3 and the lower-left components are related to the upper-right components by the 
condition = (•M|,Tree)ji- If we take the mixing to the singlet Higgs into account, the lightest 

Higgs-boson mass reads |7j : 


TO 


2 

h.NMSSM 


Mi cos 2 2/3 + X 2 V sin 2 2/3- -v 2 

Av 




sin 2/3 


+ ^ m h,l loop 


(4.5) 


As can be seen from this equation, the mixing to the singlet Higgs makes the tree-level lightest Higgs- 
boson mass smaller. The A dependence of the lightest Higgs-boson mass mainly comes from the second 
and third terms, and too large value of A makes the Higgs boson very small. There are two ways to 
decrease the mixing with the singlet: One way is to assume a small A (< 0.1), and the other is to tune 
the parameters to satisfy the relatiord, 


A- 



sin2/3 = 0 . 


(4.6) 


Conditions from positive CP-even and CP-odd Higgs boson mass-squared 

The (3,3) element in the CP-odd Higgs-boson mass matrix is given as, 

(A4p) 33 = 4A Kv u v d + \A\^-XL - 3 ksA k , (4.7) 

s 

where, in the sample parameter region we study in this paper, the third term on the right-hand side gives 
the dominant contribution. Therefore, in order for the CP-odd Higgs mass-squared to be positive, we 
must have the condition, 


ksA k < 0 , (4.8) 

in the approximation that the first and second terms in Eq. (14.71) are negligible compared to the third 
term. 

Another condition is that the (3, 3) element of the CP-even Higgs-boson mass-squared matrix 

(A4|) 33 = XA \—-- f ns (A k + 4ks) , (4-9) 

should be positive: 

—4 (ks ) 2 < ksA k , (4-10) 

where we have worked in the approximation s v u ,Vd- This condition comes from the requirement that 
the singlet Higgs-boson mass-squared must be positive in the approximation that the mixing between the 
singlet and any of the MSSM Higgs doublets is neglected. Summing up, the condition which A K should 
satisfy is 


—4(ks) 2 < ksA k < 0 . (4.11) 

In the numerical analysis presented in this paper, we give A K as an input parameter at the SUSY scale. 

3 The sum of the first three terms on the right-hand side of Eq. J4.5D is not an exact expression for the lightest Higgs- 
boson mass at tree level, but just an approximation for it. Therefore, even after making the third term vanish by imposing 
the condition Eq. (14.61) . the sum of the first three terms does not necessarily agree with the exactly maximal value of the 
tree-level lightest Higgs-boson mass obtained by diagonalizing the full 3x3 Higgs boson mass matrix, Eq. a, but only 
approximately. However, when the two conditions i-i e ftB e fi uj, v? and 4 k, 2 s 2 ^ 2A/x e fju — (B e ff + «s)usin2/3 are satisfied, 
this approximation holds with a very good accuracy. 





Constraint from non-vanishing VEV of S 


There is a condition on the model parameters from the requirement that the singlet Higgs S has a non-zero 
VEV, (S) = s ^ 0. When s v u ,Vd, the potential for S reads: 

V{S) ~ m 2 s S 2 + \kA k S 3 + ac 2 S' 4 . (4.12) 

o 

If we require that this potential has a minimum at S' = s ^ 0, and that the value of V(S) at S = s is 
smaller than V(0), we obtain the condition [7], 

Al > 9 m 2 s . (4.13) 


Constraint from Perturbativity of A 

The tree-level Higgs boson mass becomes larger for larger value of A unless we take the mixing with the 
singlet into account. However, there is a limit on the size of A which comes from theoretical consideration. 
Namely, in order for A not to blow up below the GUT scale, the value of A at the SUSY scale must be 
smaller than ~ 0.7 [7]. 

Condition from the SM-like lightest Higgs boson 

In this paper, we identify the lightest CP-even Higgs boson as the Higgs boson discovered at the LHC pQ. 
The properties of the discovered particle such as the decay branching ratios are known to be consistent 
with those of the Higgs boson in the minimal SM. This means that we have to require that the lightest 
CP-even Higgs boson in the model we consider should not be singlet-like but like the lightest Higgs boson 
in the MSSM which is known to become SM-like in the decoupling limit. 

5 Numerical Results 

In this section, we give our numerical results. 

First, we explain how we choose independent input parameters. To maximally keep the similarity 
to the cMSSM, we choose tan/3 at the SUSY scale and mo, M \/ 2 and A 0 at the GUT scale as input 
parameters. In addition, since the parameter A directly enters in the expression for the lightest Higgs- 
boson mass, we choose A at the SUSY scale as input. If we further choose either k or A\ as input, we 
can use the two tadpole conditions Eqs. mm and (14.21) to determine /i e s(= As) and H e ff(= ns + A\), 
and then use Eq. (14.31) to fix m 2 s by using the value of A K as an additional input. Below we consider 
two cases: in one case we choose k at the SUSY scale as input, and in the other case we take A\ at the 
GUT scale as input. Summing up, we consider two sets of input parameters. In one case, we choose the 
parameters below as input, 


tan/3 , X , k , A k at the SUSY scale , 

mo , Mi /2 j A 0 at the GUT scale , (5.1) 

which we call the case 1, and in the other case, we take the parameters below as input: 

tan/3 , A , A K at the SUSY scale , 

mo , Mi /2 , A 0 , A\ at the GUT scale , (5.2) 

which we call the case 2. 


Case 1 

In this case, we determine the parameters s = fi e s /A and A\ = B e s — ns by using the tadpole conditions. 
If we are to use Eq. (E2D), we have to tune k to satisfy Eq. (1TO . The value of k in this case is 


a n A 

sin 2/3 2s 


(5.3) 
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(a) k = 0.09 (b) k = 0.05 

Figure 1: Our numerical results on Br(/x —> eq) and the lightest Higgs boson mass in the MSSM-like semi- 
constrained NMSSM, case 1. The diagonal gray contours are the contours for the predicted branching 
ratios of /x — > eq. The diagonal solid and dotted red lines are the current limit and the near-future 
expected sensitivity of Br(/x —> eq), respectively. The current limits of the /x —>■ 3e and /x-e conversion 
rates are also shown by the corresponding values of Br(/x —t eq) as the solid light-blue and dark-blue 
lines, respectively. The near-future expected reach of Br(/x —> 3e) is shown by the dotted light-blue lines, 
and that of the /x-e conversion at COMET and Mu2e is shown by the dotted dark-blue lines. The gray 
line is the contour for the 126 GeV Higgs boson mass. The regions between the two red curves, the two 
solid green curves and the two dotted green curves are the areas where the Higgs boson mass is in the 
ranges, 125-127, 124-128 and 120-130 GeV, respectively. We assume mo = ATl /2 in the figures. The input 
parameters at the SUSY scale are A = 0.1 ,A K = —50 (GeV) and we take A 0 = —500 (GeV) at the GUT 
scale. The right-handed neutrino Majorana masses are taken to be M v = 5.0 x 10 14 (GeV). 


This equation means that for large tan /3 and for large A, the k parameter becomes too large, and then A 
at the scale higher than the weak scale becomes too large to be perturbative, and eventually it blows up 
below the GUT scaled We therefore do NOT assume Eq. (14.61) for the case 1, and assume small A (~ 0.1) 
to make the mixing of the MSSM Higgses with the singlet Higgs smaller, in order not to decrease the 
tree-level Higgs-boson mass. 


Numerical Results 


Our numerical results for Br(/x — > eq) and the Higgs boson mass in the case 1 are given in Figs. |T] (a) 
and (b). In the figures (a) and (b), n at the SUSY scale is taken to be 0.09 and 0.05, respectively. The 
rest of the input parameters are taken to be the same in the two figures, and the input SUSY parameters 
are A = 0.1, A K = —50 (GeV) at the SUSY scale and A 0 = —500 (GeV) at the GUT scale. We take 
mo = Mi/ 2 , and the right-handed neutrino Majorana mass is taken to be M v = 5.0 x 10 14 (GeV). The 
reason for the choice of this value of M u is that for M„ = 5.0 x 10 14 (GeV), the largest neutrino Yukawa 
coupling becomes 0(1), as can be seen from Eq. (12.231) . For smaller values of M„, the LFV rates become 
smaller since the off-diagonal entries of the slepton mass matrix become smaller, see Eq. (13.131) . 

In the figures, we plot contours for constant values of Br(/x —> eq). Since Br(/x —> eq) is roughly 

4 For A = 0.3 and tan/3 = 3, K is approximately k ~ 0.5 (the second term of Eq. I I 5.3 II is small (typically 0(0.1) or less) 
for large part of our sample parameters). The RGEs for A and k are 


nd 

167T 2 —A = A 
dt 

o d 

1671 - K, = K, 

dt 


2|«| 2 + 4|A| 2 + 3Tt(Y^Yu) + 3Tr(Y 4 Y D ) + Tr (y£y e ) + Tr(Y^. Yjy) - 3 g\ - 
6|k| 2 + 6|A| 2 . 



(5.4) 

(5.5) 


If we assume iTTeli . then a large A induces a large k, via RGEs, and A can develop the Landau pole below the GUT scale 
depending on the parameters. For small tan/3, the large top Yukawa coupling makes the right-hand side of the RGE for A 
large, and this makes it easier for the Landau pole for A to occur. 
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Figure 2: The numerical results for Br(/x —> ey) and the lightest Higgs boson mass in the MSSM +vr 
model. The diagonal gray contours are the contours for the predicted branching ratios of /j, —> ey. The 
meanings of the diagonal red/liglit-blue/dark-blue solid/dotted lines, gray line, red curves, two solid green 
curves, and two dotted green curves are the same as in Fig. H We assume too = M 1 / 2 in the figure. 
We take Aq = —500 (GeV) at the GUT scale and assume /i > 0. The right-handed neutrino Majorana 
masses are taken to be M„ = 5.0 x 10 14 (GeV). 


proportional to tan 2 /3 /toq for Aq <C too and tan/? 1, the dependence of the contours on tan/? and uiq 
are simple. In the case of the MSSM +vr model, similar results are known in literature 

In the figures, we also show the current limit and the near-future expected sensitivity of Br (/jl —> ey) by 
the solid and dotted diagonal red lines, respectively. The current limits of the /i —> 3e and fi-e conversion 
rates are also shown by the corresponding values of Br(/r —> ey) in the figures as the solid light-blue and 
dark-blue lines, respectively. Similarly, the near-future expected reach of Br(/r —» 3e) is shown by the 
dotted light-blue lines, and that of the /i-e conversion at the COMET and Mu2e experiments is shown by 
the dotted dark-blue lines. Once the PRISM/PRIME experiment is realized, it is expected to go beyond 
the COMET/Mu2e sensitivity by about two orders of magnitude, and the full region in the figures will 
be covered for this particular choice of the input parameters. 

Also shown in Figs.|T|(a) and (b) are the contours for the lightest Higgs boson mass. From the figures, 
we find that smaller re makes the Higgs boson mass smaller. We have numerically confirmed that the 
difference in the Higgs boson mass mainly comes from the values of re, and the difference in the values of 
the other parameters like A\ are not very important for the difference in the predictions for the Higgs 
boson mass. This dependence of the Higgs boson mass on re can be understood from Eq. (ESI). Namely, 
large re makes the (3,3) element of M. 2 S Xree larger and the mixing between the MSSM Higgses and the 
singlet Higgs, which makes a negative contribution to the lightest Higgs boson mass, smaller. 

From the figures, we find that there is a parameter region which is favored from the Higgs boson mass 
measurement where the predicted value of Br(^t —> ey) is within reach of the near-future experiment even 
if too is as large as ~ 4 TeV. In addition, the near-future experiments Mu3e, COMET and Mu2e can 
probe the SUSY mass scale up to ~ 5 TeV for our sample parameters. The reach will be extended further 
if the PRISM/PRIME experiment is carried out. 

We here comment on the dependence of the Higgs boson mass on re. In the figures, we take re only 
down to 0.05. For smaller values of re, for example, re < 0.03 for A = 0.1, the Higgs boson mass sharply 
decreases for decreasing re. This sharp re dependence comes from the factor (A/re) 2 in the third term of 
the right-hand side of Eq. (14.51) . If we take smaller value of A, this sharp decrease of re happens at smaller 
value of re, and hence we can take smaller re as well. 

The difference between the above results and the result in the case of the MSSM +vr model can 
become clearer if we compare Figs, [j] (a) and (b) with the prediction in the MSSM +vr model for similar 
input parameters. In Fig. [2] we show the prediction for Br(/i —» ey) in the MSSM +vr model with the 
boundary conditions at the GUT scale, Eqs. El El together with the conditions at the GUT scale 
Mi /2 = wo and [in 2 N )ij = TOg/iy, where ( m 2 N )ij ( i,j = 1,... ,3) is the soft SUSY breaking mass-squared 
matrix for the right-handed neutrinos. We take Aq = —500 (GeV) at the GUT scale and assume /r > 0. 
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(a) A a (M G ut) — —5000 (GeV) (b) ^4 a(^gut) — —2500 (GeV) 


Figure 3: Our numerical results on Br(/x —> ey) and the lightest Higgs boson mass in the MSSM-like semi- 
constrained NMSSM, case 2. The diagonal gray contours are the contours for the predicted branching 
ratios of /x —» ey. The meanings of the diagonal red/light-blue/dark-blue solid/dotted lines, gray line, 
red curves, two solid green curves, and two dotted green curves are the same as in Fig. |I] We assume 
too = Mi / 2 in the figures. The input parameters at the SUSY scale are A = 0.1, A K = —50 (GeV) and 
we take Ho = —500 (GeV) at the GUT scale. The right-handed neutrino Majorana masses are taken to 
be M v = 5.0 x 10 14 (GeV). 


We also take the same neutrino masses and mixing parameters as in Figs. |T] (a) and (b), as well as 
the same right-handed neutrino Majorana masses which are taken to be M v = 5.0 x 10 14 (GeV). By 
comparing Figs.|T](a), (b) with Fig. O we see little difference for the prediction for Br(/x —> ey) for given 
values of tan/3 and mo (= M 1 / 2 ), but the region favored from the Higgs boson mass slightly changes. For 
example, for a fixed value of tan/3 = 10, the value of Br(/x —> ey) predicted from the fixed value of the 
Higgs boson mass at to h = 126 GeV in the MSSM +vr model is ~ 10 _1 ° for this particular choice of 
input parameters, while in the case of Figs.|T|(a) and (b), the corresponding values are ~6x 10~ n and 
~ 2 x 10 -11 , respectively. We therefore conclude that in the NMSSM +vr model, the predicted value of 
Br(/x —> ey) favored from the Higgs boson mass can slightly change from the MSSM +ur model. 


Case 2 


In this case, if we are to use Eq. flop , the value of n is determined to be, similarly to the case 1, 

a n A 

sin 2/3 2s 


(5.6) 


Similarly to the reasoning in the case 1, the equation above implies that if tan/3 or A is too large, the n 
parameter at higher scale blows up and becomes non-perturbative. Therefore, if we are to use Eq. (HU), 
we need small A and small tan /3, but this choice makes the Higgs boson mass very similar to the MSSM 
case and hence is not very interesting. We therefore do NOT use Eq. dm in the case 2, either. 


Numerical Results 

In Figs. [3] (a) and (b), we give our numerical results for Br(/x —» ey) and the lightest Higgs boson mass in 
the case 2. In the figures (a) and (b), A\ at the GUT scale is taken to be —5000 GeV and —2500 GeV, 
respectively. The rest of the input parameters are taken to be the same in the two figures, and the input 
SUSY parameters are A = 0.1, A K = —50 (GeV) at the SUSY scale, and A 0 = —500 (GeV) at the GUT 
scale. We take mo = M 1 / 2 , and M v = 5.0 x 10 14 (GeV). 

In the figures, we plot contours for the constant values of Br(/x —> ey). The behaviors are very similar 
to the case 1. Also shown are the favored regions from the Higgs boson mass and the current and near¬ 
future expected sensitivities of the /x —> ey, /x —> 3e and /x-e conversion rates, similarly to Figs. [T| (a) and 
(b). 
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From the figures, we find that only in Fig. [3] (b), there is an extra Higgs-mass favored parameter 
space at the region where tan /3 and mo (= Mi/ 2 ) are both large. This difference between the two figures 
mainly comes from the difference in the value of re, and the differences in the other parameters like A\ 
enter only indirectly through the value of re in the prediction for the Higgs boson mass. 

We now explain why the changes in the input value of A\ at the GUT scale affect the value of re at 
the SUSY scale. 

To do so, we first explain the dependence of re on mo (= M 1 / 2 ) and tan/3 with fixed value of A a (Mgut)- 
Below we will show that re becomes smaller for larger mo and for larger tan /? at the region tan /3 1 in 

the parameter space shown in Figs. [3] (a) and (b). At the upper-right region of Fig. [3] (b), the value of re 
becomes re < 0.03, where the Higgs boson mass decreases relatively quickly for decreasing re, as discussed 
at the end of the discussion for the case 1 in this sectioijf]. In the parameter region shown in Fig. [3] 
(a), the value of re is larger than 0.03, and hence this relatively fast decrease does not happen. Then we 
have to explain why re is smaller in Fig. [3](b). This is because A a (Msusy) is larger in Fig. [3] (b) since 
Aa(Mgut) is larger. The relation between A a (Msusy) and re are given by re = (B e g — A\)/s and hence 
larger A\ means smaller re. 

Let us discuss the change in re for different mo (= Mi/ 2 ) for fixed values of tan/3 and A a (Mgut)- 
The value of A\ at the SUSY scale is given by solving the RGE, 

16^ 2 ^A a = 4|re| 2 A K + 8|A| 2 A a + 6Tr (Y^Tu) + 6Tr (Y^T D ) + 2Tt(Y^T e ) + 2Tr(Y^T JV ) 

+ 6^M 2 + \g\Mx . (5.7) 

5 

For our sample parameters, A a (Msusy) becomes largei@ for larger mo (= Mi/ 2 ) and fixed tan/3. There¬ 
fore, for a fixed value of tan /3, larger mo (= Mi / 2 ) makes re smaller through the relation, re = (B e g —A a ) /s. 

Next, we discuss the dependence of re on tan/3, fixing the values of mo (= Mi/ 2 ) and A\ at the GUT 
scale. Since here we are mainly interested in the difference at large tan /3 region, in this paragraph we 
assume tan/3 1. For large tan/3, A a (Msusy) becomes larger for larger tan/3 since the fourth term of 
the right-hand side of Eq. UEZD, which involves the bottom Yukawa coupling, becomes more important. 
This increase in A a (Msusy) for larger tan/3 makes re smaller for fixed mo since re = — A\)/s. 

Another reason which makes re smaller for larger tan/3 comes from the values of and B e g, although 
this effect is less important for large tan/3. The values of fx e g and B e g at the SUSY scale are obtained 
by solving the tadpole conditions, and the solutions at the tree-level are, 


Meff = 


B e ff = 


fn 2 Hd — m 2 H tan 2 /3 


tan 2 /3 — 1 


- ~ 2 Ml , 




(m T m Hd -(- v A ) sin 2/3 . 


(5.8) 

(5.9) 


Both /r e ff and B e g become smaller for larger tan/3 for our sample parameters. From the relation fi e g = As, 
a smaller fi e g means a smaller s for fixed A. From re = (B e g — A\)/s , the variation of re comes from that 
of s (= n e g/ A) and that of B e g. For our sample parameters, since the decrease in B e g due to increase in 
tan/3 has a larger effect on re than that of s, re becomes smaller for larger tan/3. 

As for Br(/r —> ey), also in the case 2, we find that there is a parameter region which is favored from 
the Higgs boson mass and in which the predicted value of Br(/i —> ey) is within reach of near-future 
experiment even if mo ~ 4 (TeV), which has not yet been probed at the LHC. 


6 Summary and Discussions 

In this paper, we have studied the cLFV in the semi-constrained NMSSM+z/r model, taking into account 
the recent results on the Higgs boson mass determination. We have considered the boundary conditions 

5 If k is too small, the (3,3) element of the Higgs boson mass matrix becomes very small and the lightest Higgs boson 
becomes singlet-like. In all the parameter regions we consider in this paper, the lightest Higgs boson is MSSM-like. 

6 At first glance, it appears that larger mo (= Mi/ 2 ) makes the gaugino masses larger, which makes the right-hand side 
of Eq. (15.711 becomes positive, and that A\ at the SUSY scale becomes smaller. However, in reality, the contributions from 
At and A 5 makes negative contributions, and the balance between the gaugino mass contributions and the At and A 5 
contributions determines the scale dependence of A\. 
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at the GUT scale to be MSSM-like and semi-constrained in the sense that the SUSY breaking parameters 
A\,A K ,rrig which are specific to the NMSSM are not necessarily equal to Ao,A 0 ,TOq, respectively. We 
have considered two cases: in one case the parameters ( s , A \, m 2 s ) are determined from the tadpole 
conditions, which we call the case 1, while in the other case (s, n, m|) are determined from other input 
parameters, which we call the case 2. 

One of the advantages of the NMSSM is that the tree-level lightest Higgs boson mass can be taken to 
be larger than that of the MSSM by taking a large value of A. In addition to this effect, there is another 
new effect in the Higgs boson sector of the NMSSM, namely, we also have to take into account the mixing 
with the singlet Higgs. This mixing can decrease the Higgs boson mass depending on the parameters. In 
the semi-constrained scenario we have considered, we find it is difficult to realize both large A and small 
mixing with the singlet at the same time. Hence in this paper we have assumed a small A (~ 0.1) which 
makes the mixing with the singlet small. 

In the case 1, we have obtained the results similar to those in the MSSM + vr model. We have also 
shown that the Higgs-boson-mass favored parameter region depends on the value of k. As the case 2, we 
have considered the case where the k parameter is not an input parameter but is a parameter determined 
from other parameters via the tadpole conditions, and we have obtained a partly different favored region 
from the case 1. In both cases, we have shown that in the NMSSM+j/r model there is a parameter region 
in which the predicted value of Br(/z —> ey) is so large that the n ej decay can be observable at the 
near-future experiment even if the SUSY mass scale is about 4 TeV. The reach will be improved further 
by the near-future experiments, Mu3e, COMET, Mu2e and PRISM/PRIME. 

Several comments are in order. In this paper we have taken the input SUSY mass parameter mo(= 
M 1 / 2 ) as high as > 1 TeV. This choice makes the squarks and gluino as heavy as multi-TeV, whose 
possibility is still not excluded by any experiments including the LHC [T] . The price we have to pay 
is that it is difficult to explain the muon g — 2 anomaly in terms of SUSY if we take the multi-TeV 
SUSY particle mass scenario, and that we have to introduce the so-called “little hierarchy” between the 
weak scale and the SUSY scale. In particular, some fine-tuning is necessary in order to keep the Higgs 
boson mass protected from large radiative corrections. Nevertheless, in SUSY models the cancellation 
between the bosonic and fermionic loop contributions to the Higgs-boson mass-squared is automatic at the 
scales much higher than the SUSY breaking scale, which decreases the degree of fine-tuning significantly 
compared to the non-SUSY minimal standard model and makes SUSY models still attractive. 
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